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ABSTRACT 

We explore the X-ray properties of a subset of the optically selected SDSS cluster sample of Goto et 
al. (2002), by analysing seven public XMM- Newton pointings, with exposure times ranging from ~ 4 to 
46 ksec. There are in total 17 SDSS clusters out of which only eight are detected at X-ray wavelengths 
with fo.5-2keV<^ 1-2 x 10~ 14 ergs cm~ 2 s _1 . For the remaining 9 SDSS clusters we estimate their 3cr 
luminosity upper limits (corresponding to L X H 5 x 10 42 ergs/sec in the 0.5-2 keV band). This relatively 
low luminosity suggests that if real structures, these galaxy aggregations correspond to poor groups of 
galaxies. Using the SDSS photometric catalogue we also derive the cluster optical r-band luminosities. 
The resulting scaling relations (L pt — L x , L opt — T x ) are consistent with those of other recent studies. 

Subject headings: galaxies: clusters: X-rays, general - large-scale structure of universe 



1. INTRODUCTION 

The intense recent interest on galaxy cluster studies 
is based on the tight cosmological constraints that they 
can provide, not only from their large-scale distribution, 
baryon fraction, clustering and dynamics, but lately also 
from their internal dynamics, its evolution and non-linear 
physical processes, that give rise in a variety of phenom- 
ena related to their formation processes (eg. Peebles et 
al. 1989; Borgani & Guzzo 2001; Rosati, Borgani & Nor- 
man 2002). However, in order for such studies to provide 
the tight constraints on cosmological parameters, that we 
have recently seen from the amazing CMB experiments (cf. 
Boomerang, Arceops, Wmap etc) it is essential to un- 
derstand and eliminate all possible systematic biases that 
enter in the construction and identification procedures of 
cluster samples (eg. Nichol 2002). 

To this end a large effort has recently been put forward 
to construct large, bias-free, samples of clusters of galax- 
ies spanning a wide redshift range (eg. Postman et al. 
1996; Olsen et al. 1999; Gladders & Yee 2000; Goto et 
al. 2002; Bahcall et al. 2003). A striking result emerging 
from such studies is that the identified clusters strongly 
depend not only on the detection algorithm but also on 
the wavelength used. Optical selection methods, based on 
galaxy overdensities in two-dimensions are biased due to 
projection effects, although more sophisticated methods, 
using color information and/or structural cluster features, 
provide a more uniform selection of clusters (cf. Postman 
et al 1996). However, even such algorithms may suffer 
from projection effects and biases towards more evolved 
galaxy populations. Alternatively, X-ray selected cluster 
samples are less prone to biases due to the centrally peaked 
X-ray emitting cluster core which can be seen even at large 
[z ~ 1) distances (eg. Stocke et al. 1991; Castander et 
al. 1995; Ebeling et al. 1996a, 1996b; Scharf et al. 1997; 
Ebeling et al. 2000; Bohringer et al. 2001; Gioia et al. 
2001; Rosati, Borgani & Norman 2002). Nevertheless even 
such identification procedure could be biased towards the 



most virialized (and hence more X-ray luminous) clusters. 
A novel cluster detection approach has recently been pro- 
posed using simultaneous multiwavelength data with the 
use of Virtual Observatory data (Schuecker et al. 2004). 
A generic problem, however, arises from our ignorance of 
the cluster formation and evolution details which further 
complicates the characterization of clusters found at high 
redshifts and their use as cosmological probes. 

Attempts to understand the systematic biases that en- 
ter in cluster selection procedures have been recently pre- 
sented in the literature (cf. Donahue et al 2002; Basilakos 
et al 2004 and references therein) . This letter further con- 
tributes to this issue by investigating the X-ray properties 
of a subsample of the Goto et al (2002) clusters detected 
in the Sloan Digital Sky Survey (SDSS) with the Cut 
and Enhance method (CE), which is based on galaxy col- 
ors 1 . We correlate the positions of these CE cluster can- 
didates with the extended X-ray source position, in order 
to investigate which of the clusters show significant X-ray 
emission, and thus have a larger probability of being true 
aggregations of galaxies. Throughout this paper we adopt 
the concordance cosmological model (f2 m = 1 — 17a = 0.3), 
although all relevant quantities are partametrized accord- 
ing to H — 100 h km s~ 1 Mpc _1 , in order to facilitate 
comparisons with other studies. 

2. xmm-newton data analysis 

We use XMM-Newton archival observations, with a pro- 
prietary period that expired before September 2003, that 
overlap with the first data release of the SDSS (DR1; 
Stoughton et al. 2002). Only observations that use the 
EPIC (European Photon Imaging Camera; Striider et al. 
2001; Turner et al. 2001) cameras as the prime instrument 
operated in full frame mode were employed. For fields ob- 
served more than once with the XMM-Newton we use the 
deeper of the multiple observations. A total of 7 fields are 
used with Galactic Nh in the range 2 — 8 x 10 20 cm~ 2 and 
PN good time intervals between 4 and 46 ksec (see Table 
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1). Note that two of the fields have clusters as prime tar- 
gets (Abell 267 and RXJ0256. 5+0006). The X-ray data 
reduction, source detection and flux estimation are carried 
out using methods described in Georgakakis et al. (2004) . 

In total there are 17 CE optical cluster candidates in 
the region covered by the analyzed XMM-JVewton obser- 
vations (see Table 2). We search for extended X-ray 
sources around the 17 CE clusters, in the 0.5-2 keV merged 
PN+MOS images, where available, using the emldetect 
task of SAS. We find that extended sources are associ- 
ated with 8 out of the 17 clusters (note that two of these 8 
clusters were the prime targets of two of the XMM-Newton 
pointings). These were visually inspected to exclude the 
possibility of spurious detections due to their possible as- 
sociation with the CCD gaps, hot pixels or the field of 
view edges. Two examples of poor clusters (No 2 and 16 
of Table 2) are shown in Figure 1, where the SDSS r-band 
image is overlayed with the X-ray contours. 

Note that 3 out of the 17 CE optical clusters are also 
found in the Bahcall et al. (2003) cluster candidate list, 
which is based on the joint application of two different al- 
gorithms on the SDSS photometric galaxy data, while it 
contains 2 more clusters not found in Goto et al. Out of 
the 3 common clusters, two are detected in X-ray's (No. 5 
and No. 13 of our Table 2). 

For each cluster candidate we estimate the X-ray flux 
from a region which encompasses all the extended emis- 
sion. Then we correct to the total flux by integrating a 
King profile to infinity: l x oc [1 + (r/r c ) 2 ]~ 3 ^ +1 / 2 , where 
(3 is the ratio of energy per unit mass in galaxies to that 
in gas. We fix the exponent of the profile to (3 = 0.5 and 
the core radius to 100 kpc for clusters with (uncorrected) 
luminosities L x < 10 43 h~ 2 ergs/sec and 200 kpc for those 
with L x > 10 43 h~ 2 ergs/sec. We find that varying the 
core radius by ±50 kpc changes the total flux by ~ 8 per 
cent. The conversion from count rates to flux is performed 
assuming a Raymond- Smith (RS) spectrum with the tem- 
perature derived from the spectral fit (see below) . For the 
clusters that have not been detected we have estimated 
their 3er upper limits, by extracting the counts in a 1 ar- 
cmin circle around the cluster optical center. Background 
counts were estimated from nearby regions. Then count 
rate upper limits are evaluated using the method described 
in Kraft et al. (1991). We convert the count rates to flux 
assuming a RS spectrum with a temperature of 2 keV, typ- 
ical of poor clusters of galaxies. The total flux is estimated 
using the King profile above. 

We investigate the X-ray spectral properties of the eight 
detected clusters by performing individual spectral fittings 
using the XSPEC vll.2 package. For four clusters which 
have enough photon statistics (5, 11, 16, 17) we have used 
the standard y 2 analysis. For the other four we use the 
C-statistic instead (Cash 1979) which is proper for fitting 
spectra with a limited number of counts. We fit simulta- 
neously the PN and MOS data in the energy range 0.3-8 
keV. We use a RS spectrum with the abundance fixed at 
Z = 0.3, absorbed by the Galactic hydrogen column den- 
sity listed in Table 1. 

Finally we present in Table 2 the optical and X-ray prop- 
erties for our clusters. Specifically, we list their optical 
equatorial coordinates (J2000), the detected flux (bold) or 
the 3(7 upper flux limits in units of 10~ 14 crg cm~ 2 s _1 , the 



cluster redshift (spectroscopic or photometric estimated by 
the CE method) , optical r-band luminosities (see next sec- 
tion), X-ray luminosities (erg s -1 ) in the 0.5-2 keV band, 
best-fit temperatures and x 2 values, where available, of 
the spectral fits. 

It is interesting that the majority of the CE optical clus- 
ter candidates are not identified in X-rays which could be 
either because they are poor, dynamically young, systems 
which have very low X-ray emission, or the result of pro- 
jection effects. The 3<t upper X-ray luminosity limits are 
low, for most cases L .5-2kcV < 3 x 10 42 h~ 2 erg s _1 , 
corresponding to the emission of poor groups of galaxies. 
In order to investigate this issue we determine the optical 
properties of all 17 CE cluster candidates. 

3. OPTICAL SDSS DATA ANALYSIS 

We estimate the optical luminosity of the cluster can- 
didates by cross-correlating the cluster positions with the 
SDSS photometric r-band galaxy survey. For each clus- 
ter we select all galaxies with m r < 21 (a limit to which 
there is excellent star-galaxy separation) falling within a 
distance of 0.5 h~ Y Mpc from the cluster center (for the 
cases where no spectroscopic redshift was available we used 
the estimated Goto et al. photometric redshift). How- 
ever, the cluster galaxy membership is affected by fore- 
ground/background contamination and therefore we esti- 
mate the local background by using a circular area of a 
1.5 degrees radius around each cluster center. We first 
exclude the high surface density regions and then mea- 
sure the background surface density, Ub- The background 
galaxy counts, projected on the surface covered by each 
cluster at its rest-frame, are estimated by: A^ack = ^d 2 ab, 
where d is the angular cluster radius corresponding to 0.5 
h^ 1 Mpc at the cluster rest-frame. Then the predicted 
number of the real members for each particular SDSS clus- 
ter is found by: 

N r = N ohs ~ N hac . (1) 
We assume that for each cluster the luminosity function of 
galaxies can be described by a double Schcchtcr function 
(eg. Smith et al. 1997), ie., *(L) = + 2/ 2 (L)], 

with 

f i (L) = ^(-^-) a, e X p(~) i = l,2 (2) 

where L M ~ 3.21 x 10 11 h~ 2 L with a x ~ -1.0 and 
L*, 2 ^ 2.02 x 10 10 h- 2 L with a 2 ~ -1.7, respectively. 
Thus, the number of cluster galaxy members N T can also 
be written: 

poo 

N r = V clus / <S>(L)dL , (3) 

where L m i n (r) is the minimum luminosity of a galaxy at 
distance r (corresponding to the flux limit), r{z) is the lu- 
minosity distance and V c i us is the volume occupied by the 
cluster. 

From eq.(l) and eq.(3) we can determine for each clus- 
ter the unknown value of normalization constant <f>„ and 
thus obtain the total cluster luminosity by: £ pt = 
Vdus/^ . r r \L§>(L)dL. Combining the above system of 
equations we obtain the final expression for L opt : 
g°. {) L[f 1 (L) + 2f 2 (L)}dL 

^ ~ Nl /L°l (l) [/iW + 2/ 2 (i)]dL " 5 " ; 
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The optical cluster luminosities and their 2a Poisson un- 
certainties are presented in Table 2. We find that one clus- 
ter, also undetected in X-rays, has a negative N r value and 
is therefore suspect of being a fluke. 

Under the assumption that both, cluster gas and galax- 
ies are in hydrostatic equilibrium with the underline grav- 
itational potential and also that the clusters have an ap- 
proximately constant M/L, one expects power-law scaling 
relations between L opt , L x and T x . Although our sam- 
ple is quite small we present as an example in Figure 2 
the L opt — L x relation and the best fit log-log relation 
(solid line). The dashed lines encompass the correspond- 
ing Popescu et al. (2004) relation. The cluster detections 
are shown as solid points and those with only upper limits 
in L x are shown as vectors, while we also plot their 2a 
L opt uncertainties. Note that the cluster with the high- 
est 3(7 upper L x limit (No 8 in Table 2) is located on the 
XMM- Newton pointing with the lowest exposure time (3.8 
ksec) and thus we are unable to provide tight constraints 
on its X-ray properties. 

Performing a linear regression in log-log space between 
the different observational quantities we find the follow- 
ing relations: L opt /L & = io"-«>(±o.io) L o.43(±o.i3) and 
L ™t/L = 10 10 - 8 °( ±0 - 15 ) r i-«»(±o.3i) ) with excellent fits 
(jq ~ 0.97 and 0.7 respectively). The above fits are in 
good agreement (within the la errors) with the Popescu 
et al (2004) relations, derived from a cross-correlation of 
the RASS and the SDSS (114 clusters). Furthermore, 
they are in relatively good agreement with the theoreti- 
cal expectations for virialized structures (L opt oc L x 5 and 
L opt oc T x b ). 



4. CONCLUSIONS 



In order to investigate the X-ray emission of the Goto 
et al. (2002) SDSS clusters, we have analysed seven public 
XMM-Newton pointings with exposure times ranging from 
4 to 46 ksecs. Out of the 17 such SDSS clusters, found in 
the regions covered by these XMM- Newton pointings, only 
eight (~ 47%) are found to have extended X-ray emission 
with limiting flux f x <~ 1.2 x 10~ 14 ergs cm~ 2 s _1 . Note 
that two of these eight SDSS clusters, which are also the 
most X-ray luminous ones, were the prime targets of two of 
the XMM- Newton observations analysed. The remaining 9 
clusters are probably either very poor and/or dynamically 
young clusters with weak X-ray emission or the results of 
projection effects. The derived L pt — L x and L pt — T x 
relations, using the clusters with detected X-ray emission, 
are consistent with those of other recent determinations. 

The present work provides a first glimpse of the obser- 
vations that the future Duo X-ray mission will perform. 
It will survey a wide angle (6000 sq. deg) in the 0.3-10 
keV energy band with the aim of using <~ 10000 clusters 
to map the geometry of the Universe and determine with 
high precision its Dark Energy component. Duo will tar- 
get areas covered by the SDSS in order to facilitate the 
cluster rcdshift determination. Our work helps to deter- 
mine the X-ray flux depth required to identify a significant 
fraction of the optically selected cluster population. 

This research was jointly funded by the European Union 
and the Greek Government in the framework of the 
project 'X-ray Astrophysics with ESA 's mission XMM- 
Newton '. within the program 'Promotion of Excellence 
in Technological Development and Research' and by the 
Greek-British scientific bilateral program (2002-2003) en- 
titled "Observations of clusters and groups of galaxies 
with XMM-Newton". We thank the XMM- Netwon science 
archive for the provision of the X-ray data used. 
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Fig. 2. — The L op t — L x relation from the cluster X-ray detections. The dashed lines delineate the Popescu et al. (2004) results, while the 
thick continuous line corresponds to our best fit. We also show the clusters with upper L x limits (vectors). 
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Table 1 

The archival XMM-Newton pointings used in this study. 



RA 


Dec 


Filter 


N H 


PN cxp. 


MOS1 cxp. 


Field name 


Obs. ID 


P.I. 


(J2000) 


(J2000) 




(10 2() cm" 2 ) 


(sec) 


(sec) 








00 43 20 


-00 51 15 


Medium 


2.33 


15 700 




UM 269 


0090070201 


Reeves J. 


01 52 42 


+01 00 43 


Medium 


2.80 


5 800 


17200 


ABELL 267 


0084230401 


Kneib J.P. 


01 59 50 


-00 23 41 


Medium 


2.65 


3 800 




MRK1014 


0101640201 


Aschenbach B. 


02 56 33 


-00 06 12 


Thin 


6.50 




11600 


RXJ0256. 5+0006 


0056020301 


Arnaud B. 


03 02 39 


+00 07 40 


Thin 


7.16 


38100 


46 900 


CFRS 3H 


0041170101 


Gear W. 


03 38 29 


+00 21 56 


Thin 


8.15 


8 900 


6 700 


SDSS 033829.31+00215 


0036540101 


Brandt W. 


12 45 09 


-00 27 38 


Medium 


1.73 


46 300 


55 500 


NGC 4666 


0110980201 


Jansen F. 



Table 2 

x-ray flux upper limits and detections (in bold) of the goto et al. (2002) ce cluster candidates. 
Photometric and spectroscopic redshifts are indicated by the superscript 1 and 2, respectively. The second 
column lists the corresponding number of the bahcall et al. (2003) sdss clusters. 



# 


^Bahcall 


a 


6 


fx 


z 


L/L Q (*10' L ) 




T(kcV) 




1 




10.89 


1.015 


8.3 


0.195 1 


0.9 ±0.4 


42.60 


I O+U.3 

-■-• o -0.2 




2 




10.96 


0.774 


3.7 


0.311 2 


0.9 ±0.6 


42.72 


1.418:5 




3 




10.72 


0.726 


4.3 


0.265 2 


1.9 ±0.8 


42.62 


2.9+™ 




4 


173 a 


28.05 


1.108 


1.7 


0.243 2 


1.6 ±0.7 


42.07 




5 


174 


28.18 


0.999 


174 


0.230 1 


4.0±1.0 


44.06 


6 9+ ' 3 

O.»_ 2 


997/768 


6 




28.18 


1.138 


2.7 


0.219 2 


1.3 ±0.6 


42.24 




7 




28.24 


1.055 


3.6 


0.231 2 


1.2 ±0.6 


42.40 






8 




29.84 


0.449 


2.1 


0.402 2 


0.9 ±0.7 


43.19 






9 




29.88 


0.215 


1.5 


0.367 2 


1.1 ±0.7 


42.43 






10 




29.88 


0.283 


1.8 


0.356 2 


0.3 ±0.4 


42.55 






11 




44.14 


0.106 


32.4 


0.360 1 


3.5 ±1.3 


43.83 


r a+0.7 
D - D -0.5 


202/103 


12 




45.82 


-0.005 


0.8 


0.277 2 


0.6 ±0.4 


41.90 






13 


259 


54.49 


0.485 


1.2 


0.323 1 


2.7 ±1.0 


42.28 


1 r+0.7 




14 




54.54 


0.274 


6.6 


0.186 2 


0.8 ±0.4 


42.45 




15 




191.13 


-0.527 


0.9 


0.231 2 




41.80 






16 




191.22 


-0.560 


6.9 


0.232 1 


1.1±0.5 


42.70 


1 6+ ' 2 

-■-• D -0.2 


128/125 


17 




191.22 


-0.445 


4.3 


0.231 1 


O.6±0.4 


42.50 


1 r+0.3 


168/115 



a at a distance of ^2.3 arcmin (<~ 0.5 h 1 Mpc at the cluster rest-frame). 



